Electromigration of Single-Layer Clusters 
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Single-layer atom or vacancy clusters in the presence of electromigration are studied theoretically 
assuming an isotropic medium. A variety of distinctive behaviors distinguish the response in the 
three standard limiting cases of periphery diffusion (PD), terrace diffusion (TD), and evaporation- 
condensation (EC). A general model provides power laws describing the size dependence of the 
drift velocity in these limits, consistent with established results in the case of PD. The validity 
of the widely used quasistatic limit is calculated. Atom and vacancy clusters drift in opposite 
directions in the PD limit but in the same direction otherwise. In absence of PD, linear stability 
analysis reveals a new type of morphological instability, not leading to island break-down. For 
strong electromigration, Monte Carlo simulations show that clusters then destabilize into slits, in 
contrast to splitting in the PD limit. Electromigration affects the diffusion coefficient of the cluster 
and morphological fluctuations, the latter diverging at the instability threshold. An instrinsic 
attachment-detachment bias displays the same scaling signature as PD in the drift velocity. 
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I. INTRODUCTION 

Recognized as a key source of size limitation in electri- 
cal devices, surface electromigration has important. con- 
sequences for surface morphology. Latyshev et aM dis- 
covered that it can induce step bunching on Si(lll) vic- 
inal surfaces. How electromigration affects more com- 
plex surface structures is still poorly understood. In 
this paper we explore the effect of electromigration of 
atoms or vacancies on single-layer clusters. Responses of 
these clusters to electromigration sharcpsimilarities with 
void behavior in metallic electric linesEI and with elec- 
tromigratiijai of oxygen disordered domains in YBaCuO 
thin filni^la and lies in the active area of driven diffusive 
systems Jj 

In addition to the intrinsic interest in how island prop- 
erties are affected by the symmetry-breaking electromi- 
gration force, we will point out further mesoscopic ways 
to address current controversies about the dominant 
mechanism of mass transport.B Moreover, these proper- 
ties are needed to model surfaces undergoing electromi- 
gration at larger scales. For example, electromigration- 
induccd coalescence models should use non-equilibrium 
diffusion constants and steady-state drift velocities. 

In the presence of an electric current, the electromi- 
gration force is usually described as F = z*eE, where e 
is the magnitude of the electron charge and z* is an effec- 
tive valence, generally non-integer, which takes into ac- 
count both the electrostatic interaction between the elec- 
tric field E and the charge distribution on the affected 
atoms ("direct" force) and the frictional force resulting 
from the transfer of momentum from the charge carri- 
ers to these atoms ("wind" force). The long-standing 
controversies regarding bulk electromigration highlight 
the complexity of microscopic calculations, and partic- 
ularly how conduction electrons may screen the direct 
charge of the atoms.Q Non-trivial aspects of the flow of 
the electron cloud about the diffusion path might also 



contribute to the direct effect on the surfaccH For met- 
als, the electromigration effect is primarily due to wind 
force; calculated values include z* « —30 for Al on scnu=. 
infinite jelliumEl and z* « -21 for Cu on Cu 
(Near a step edge, one can expect quantitative but not 
order-of-magnitude changes in z*; e.g. for an atom in 
a close- packed step edge on Al{001}, z* w -430) For 
semiconductors, bath forces are small, with a resultant 
\z*\ « 0.001 — 0.1 JIj a characteristic length can be as- 
sociated with electromigration: ^ = ksT/F. In typical 
experimental conditions, ^ ~ 10® atomic spacings for Si 
and f ~ 10^ for metals. We base our study of the equi- 
librium behavior of these islands on the framework of the 
continuum theory of Khare et aZ£3 We consider electro- 
migration to be a perturbation inducing a macroscopic 
current, as done previously in Refs. ^ and |lj. This ap- 
proximation is justified because ^ is much larger than 
the atomic spacing. The adatom flux on terraces is then 
equal to c{v). The mean drift velocity {v) of an adatom 
due to electromigration is calculated from the Einstein 
relation: 



{v)=D 



F 



D 
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Single-vacancy motion may dominate over atom mo- 
tion in the mass transport op .same surfaces. (There is 
evidence of this on Cu {001}.t£lE3) In that case one must 
use an effective charge z* appropriate for vacancies and 
reverse the signs of step curvatures in computing restor- 
ing forces due to line tension. The resulting modifica- 
tions are relatively straightforward but tend to muddle 
the subsequent descriptions. Accordingly, we do not ex- 
plicitly consider vacancy transport below. 

In the notation of Ref. |l^, there are three limiting 
mass transport modes: periphery diffusion (PD), ter- 
race diffusion (TD), and (two-dimensional) evaporation- 
condensation (EC) or attachment-detachment (see 
Fig. 1). In particular, there is a distinctive size depen- 
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dence of the tracer diffusion constant D^'^ of a cluster 
at thermodynamic equiUbrium, with characteristic ex- 
ponent a =3, 2, and 1, respectively: 



n't 



(2) 



where Rq is its average radius. Intermediate values can 
also be obtained in restricted regions of parameter space 
as one crosses over from one limiting regime to another. 
As a first qualitative approach, the mean cluster drift 
velocity should be given by the Einstein relation in a 
way very similar to the drift velocity of one adatom (Eq. 
(l|)): it is equal to the cluster mobility (or equilibrium 
diffusion constant I?^*) time the force exerted on the 
whole cluster F = irR^F/a? , and divided by fcsT: 



\V\ « Dl'^ 



F 



Rn 



(3) 



This scaling law will be confirmed later. The interme- 
diate expressions are in fact exact for models in which 
atoms approach and leave the island only from the lower 
edge (one-sided step), or in the PD regime. (In the more 
general case, some subtleties arise.) 

We further investigate the shape changes associated 
with steady states. Circular steady states are found in 
the PD and TD regimes. But when attachment and de- 
tachment of atoms to the steps is not instantaneous, the 
cluster elongates. This can be understood by analogy 
to the Bernoulli effect. Stability of these steady states 
is studied. Besides the expected splitting of clusters for 
PD, we find another morphological instability in the TD 
case. Monte Carlo simulations indicate that it leads to 
slit formation. Fluctuations and the non-equilibrium dif- 
fusion constant are also calculated within the framework 
of a Langevin model, and appear to be corrected by a 
bias term proportional to 1 /^^ for weak electromigration. 



II. MODEL 

Our formal approach is inspired by a previous exam- 
ination qL non-equilibrium step meandering on vicinal 
surfaces .im For simplicity and to avoid distracting com- 
plications, we take all step and terrace properties to be 
isotropic. We forbid adsorption onto or desorption from 
the terrace. Evolution of the mobile-atom concentration 
on terraces is then given by the following conservation 
law: 



(9tc= -V- J, 



D 



J = -L'VcH c X 



(4) 
(5) 



J is the mass flux on terraces, and D is likewise the 
adatom diffusion constant on terraces. The correlations 
of the conserved noise q are given below. In the widely 
used quasistatic limit, one considers that the adatom 



concentration reaches a steady state on time scales much 
smaller than that of step motion. Then 



V- J = 



(6) 



has to be solved on terraces instead of Eq.(Q). Adatom 
electromigration is taken to be unaffected by the pres- 
ence of steps; thus, ^ is uniform. Atom exchange between 
the step and the 2D "gas" on the terrace is proportional 
to the deviation from equilibrium, expressed in the linear 
kinetic relation: 



(7) 



In our notation, -I- or — denote the lower or the upper 
side, respectively, of the step forming the boundary of 
the island. Thus, for atom islands, + is the exterior while 
for vacancy islands it is the interior. The unit vector fi 
normal to the step points toward the -I- side. The kinetic 
coefficients i'± describe attachment and detachment at 
the two sides of the step edge, and ri± are non-conserved 
noises. Attachment lengths are defined by d± = D/i'±; 
the larger d±, the smaller the chance that an atom will 
detach from or attach to the step. When d± are small, 
the dynamics are diffusion limited (TD). For large but 
finite d± , atoms attach only after a large number of tri- 
als, and dynamics are thus limited by attachment and 
detachment (EC). In the limit that d_ — > oo, the model 
reduces to an exterior (or interior) model for atom (or 
vacancy) islands, because atom exchange occurs on only 
one (viz., the lower) side of the step. The equilibrium 
concentration Ceq must include corrections due to bound- 
ary curvature, as given by the Gibbs-Thomson relation 



--eq 



= c" exp(rK) , 



(8) 



where c^^ is the equilibrium concentration in the vinicity 
of a traight step, k is the step curvature (counted positive 
for a convex step cf. Eq. (p^)), and F is the capillary 
length: 



F = a^p/ksT. 



(9) 



The step stiffness /? is taken to be isotropic, as noted 
above; a is the lattice constant, and T is the temperature. 
Finally, the normal velocity Vn (i.e., along n) of the step 
is given by mass conservation 



Vn= a n ■ (J , - J| ) - dsJs, 



(10) 



where s is the arclength along the step. We consider 
the low-concentration limit, where c ^ i.e. the 

atom concentration is taken to be much lower than in the 
bulk. Hence, the second term, representing advection, in 
J|± — —D\/c± + hVnC± is neglected. Jst is the mass flux 
along the step. In addition to the usual relaxation term 
related to chemical potential gradient, electromigration 
induces a new term D^Fstny/kBT = DLUy/^^st, with 
^st = ksT/Fst- Fst is the force exerted on mobile atoms 
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at the step edge, and Uy is the y component of n. Thus, 
we have 

J St = a^ds[R{e) cosie)] - aDLdsiVn) - q^u (11) 

where qst is a Langevin force. R{0) is the distance 
from the center of mass of the cluster, and 6 is the po- 
lar angle with respect to the x axis (see Fig.l). The 
macroscopic step diffusion constant could be defined as 
Dl = a?CstDst, where Cst is concentration of mobile 
edge-atoms, and Dst is the diffusion constant for mo- 
tion of these atoms along the step, or using the Kubo 
formula.Ea In the following, we will take Dl to be uni- 
form along the step and constant in time. 

In a local thermodynamical equilibrium approxima- 
tion, noise correlations are written 

{qi(x,y,t)qj{x',y',t')) = 2Dc{x,y,t)S , 
2c{s,t)\± 



{V±{s,t)v±is',t')) = 
{qst{s,t)qst{s\t')) 



a^Dr 



'-5{s~ s')5{t~t'). 



{dtv. 



CM) 



and the non-equilibrium tracer diffusion constant: 



|rcM- Vtp) 
it 



(15) 



(16) 



which reduces to the usual equilibrium cluster diffusion 
constant Z?^* when V = 0. This quantity could also 
be defined as half the diffusion constant describing the 
relative motion of two identical islands, whether they 
are drifting or loat. This latter definition has been useful 
experimentallyE3 Shape fluctuations of the cluster are 
measured by the island roughness 



At equilibrium, using Eq. (H), we find: 



2 _ SksTRo 

eq 
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(17) 



(18) 



As a static property, W"^^ does not depend on the trans 



.^(s _ s')5{t — t') , (12) port mechanism, just like the static spectrum in Eq.(14) 



where 5 = S{x — x')5{y — y')5{t — t')Sij, and c{s,t)\± is 
the adatom concentration on the terraces in the vicinity 
of the island edge. 

In describing the fluctuations of these islands, we con- 
sider only the case in which overhangs can be ignored. 
Thus, the periphery can be described in polar coordi- 
nates by R{d). For small fluctuations from the circular, 
R{0) = Rq + p{0). The Fourier transform of p is given 

by 



de 



dt 



2vr./-oo 27r 



p{0,t)e' 



(13) 



The equilibrium properties of such clustcES- (for _F=0) 
were calculated in detail by Khare et a/£3 They are 
able to evaluate the equilibrium diffusion constant un- 
der general conditiops. Equipartition of energy provides 
the static spectrum£j 



(|P"|^)oq 



ksTRo 



27r(n2 - l)/3 



(14) 



The divergence of (|pi|^)oq reflects the absence of an en- 
ergy cost to translate the whole island: this is a "Gold- 
stone mode." The mode n = is forbidden due to 
area conservation for PD dynamics. But generally in 
the presence of non-conserved dynamics, the system will 
tend to minimize step length, and clusters are unsta- 
ble: they expand or shrink, depending on the precise 
kinetics and environment. As an exarnple, an assembly 
of clusters will coaxsen due to Ostwalc£3 or coalescence 
( "Smoluchowski" )EiJ ripening. 

Let rcM denote the position of the center of mass of 
the cluster. We define the mean velocity: 




FIG. 1. Three processes are involved in mass transport: 
Diffusion along the step, diffusion across the terraces, and 
attachment-detachment of adatoms at the step. When one 
of these processes is slower than the others, we get the PD, 
TD, or EC regime, respectively. Electromigration force F is 
taken to be along the x axis in the text. R is the distance 
from the center of the island, 6 is the polar angle with respect 
to the X axis, and s is the arclength along the step. _Ro is the 
mean radius. 

A major goal of our analysis is to explore proper- 
ties of driven clusters on the dominant mechanism of 
mass transport. A crude criterion for being in the PD 
regime follows from checking wether the diffusion con- 
stant is dominated by line diffusion or diffusion across 
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the terrace. Comparing expressions (g5|) and ( p4[ ) for 
D^"^ given below, one finds: Rq {Rq + d)£c, where 
d = (i+ or (i_ and ic = D^/ Dac^^^K^ The opposite case, 
i?o 3> (i?o + rf)^c, will be denoted as the 2D transport 
regime. In this regime, TD and EC limits are defined 
by i?o d+ or and i?o *C and respectively. 
Note that the cluster radius i?o is involved in these re- 
lations: it is a geometrical cut-off for long-wavelength 
fluctuations. A more general definition of these regimes 
involving wavelength dependence can be found in Ref. 



III. PD LIMIT 

A. Steady states and stability analysis 

Some known results about 2D voids in metal electric 
line can be transposed directly to the case of vacancy is- 
lands in the PD limit. Best known is "Ho's law," which 
states that the drift velocity varies inversely with the 
radius.Ba Behavior consistent with this relationship was 
seen in recent Monte Carlo simulations of electromigr|ai 
tion on Cu{001} using semi-empirical energy barriers.tj 
In our formalism, the drift velocity (in the a;-direction) 
of a 2D cluster with mean radius i?o is: 



istRo 



(19) 



(j> = 1 (—1) denotes atom (vacancy) islands. Thus, atom 
islands and vacancy clusters drift in opposite directions. 
With current densities about 10^-10^ A cm~^ at 600 
K on metals, we find ^st ~ lO^A. Assuming next an 
activation barrier of ~0.5eV along the periphery and a 
concentration of mobile atoms per site c^t ~ 0.1, we 
estimate Dl ~ lO^A^/s; thus, for an island of radius 
lO^A, V « lOA/s. 




Vacancy Island 




Atom Island 





EC-TP 



FIG. 2. Drift direction 
mass transport mechanism. 
Ehrlich-Schwoebel effect. 



of island as a function of 
We have assumed a direct 



Stability analysis for voidfl'0 can also be applied to 
our case. Since we assume isotropic step properties, the 
equilibrium shape of the island is circular. With in- 
creasing electromigration, this steady-state shape per- 
sists so long as the island is stable. A hand-waving ar- 
gument allows us to retrieve the instability threshold 
found in Refs. |2^ and We simply have to com- 

pare fluxes along the step close to a protrusion of the 
step inside the cluster. The flux contributing to desta- 
bilization is Jdestab ~ I^l/S,- Stabilizing effects related 
to hue tension effect result in a flux Jstab ~ DlT/Rq, 
where curvature changes are taken to be ~ l/i?o be- 
tween points separated by a distance ^ Rq. The protru- 
sion increases if Jdestab > Jstab- The instability criterion 
is then Ixl > Xc, where x = i?o/r^ = (Pom/ l3){RoUf. 
According to a linear stability analysis, Xc = 10.65.E3'EZl 
With the preceding parameters and /3 =0.3eV/a, we find 
that islands with radius larger than i?c ~ 5 • lO'^A will 
be unstable. The instability appears at a characteristic 
time T = Xc^^^ /{^Dl) « lO^s. Thus, this phenomenon 
should be observable experimentally. The instability cri- 
terion is the same for atom and vacancy islands. 

Since electromigration of single-layer clusters is unaf- 
fected by their shape, our voids can be characterized as 
conducting. In contrast, voids in electric lines are essen- 
tially insulating. The resulting current-crowding effects 
qualitatively change the nature of the instability: the 
void is linearly stable, but becomes unstable under finite 
pertudiations. At late stages of this instability, the void 
splits.y 



B. Diffusion constant and fluctuations 



Keeping only the line diffusion term in Eq. (|10[), one 
gets in the frame moving with the mean velocity V of 
the island: 



Vn - Vn, 



-ds[Jst] 



— DLd,[Ri9)cosie)] - aDLds(TK) - q^t 
4st 



, (20) 



where is the x component of the unit vector n normal 
to the step. The normal step velocity is written as: 



dtR 



[1 + (deR/R) 



211/2 



(21) 



and the step curvature is defined with the sign conven- 
tion: 



K — (j)- 



R- 



dfjoR 



[1 + (deR/Rn 



3/2 



(22) 



We then linearize Eq. (|20|) for small deformation of the 
island p{9). With the Fourier transform of p{0,t) as 
defined by Eq. (|l|), Eq. (M) takes the form: 
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[iuT + n^(n^ - 1)] pnui 

ITIT 

-n[{n + 2)pn+i,c^ + (n - 2)pn^i,^] 4>x = -^Hnuj-, (23) 

where r = uDlT / Rq. This equation imphes that po is a 
constant; but this constant must vanish (po = 0) since p 
is defined as the departure from the mean radius Rq. In 
Eq. (e3) the eigenvalues ioj do not depend on the sign 
of X (as expected physically). Thus, these eigenvalues 
also do not depend on 0. Hence, as stated above, atom 
and vacancy islands have the same instability criterion 
Ixl > Xc- 

We first consider the case x ^ 1- We then expand step 
fluctuations in the form p„ = pi"-* + xPn'' + X^Pn^ + ... . 
To 0th order, we simply get Eq. ( p3[ ) without the term 
proportional to x- The diffusion constant is calculated 
using the Fourier space relation: 



^c = (|pi|')A. 

The equilibrium diffusion constant reads: 



nRl 



(24) 



(25) 



To higher orders in x, Eq- ( P^ ) shows the Langevin force 
q does not explicitly intervene, ultimately because Dl is 
supposed to be independent of electromigration; hence, 



= n[(n + 2)pi"^\_^ + (n- 



(m) 



0X 



(26) 



Using this relation, we calculate the first correction to 
the cluster diffusion constant: 



\ + 0{x') 



(27) 



The non-equilibrium cluster roughness can be calculated 
in a similar way. We find: 



W% (1 + bx^ + 0{x^)) 



(28) 



where h ~ 0.23 is a numerical constant resulting from 
infinite summations. Since x is of order unity, this cor- 
rection can be non-negligible. Note that the first correc- 
tion to these quantities is proportional to the electromi- 
gration force squared. Correspondingly, W and Dc are 
invariant under the inversion symmetry F —F. 

Fluctuations increase with F. In a linear theory, they 
should diverge at the instability threshold. If we define 
e = Xc — then close to the threshold, Eq. ( |2^ ) can be 
written in terms of the vector p = {p„ : \n\ > 2} 



liuop = M.p + 6, 



(29) 



where & is a noise term, I is the identity matrix and the 
matrix M is expanded as M = Mc + eMi, where Mc 



has one null eigenvalue (as can be seen, for example, in 
Ref. table II). It is then easily shown that the squared 
roughness 



= {fp) 
dco 



{b^{iujI + M)-'-''\^iujI + M)-'-b) (30) 



2tt 



diverges like e ^ when e 
instability threshold. 



0. Hence, approaching the 



.-1/2 



(31) 



Note that Dc does not diverge: the instability concerns 
the morphology of the island, not its motion. The diver- 
gence of W is analogous to the case of an isolated st^ 
subject to a morphological instability during growth,c3 
where the same exponent was found. Close to threshold, 
fluctuations become large, and nonlinear effects should 
be addressed, as shown in Ref. ^ during growth. 



C. Monte Carlo simulations 

We here use kinetic Monte Carlo simulations to show 
that the morphological instability of an island under pe- 
riphery diffusion leads to splitting. 

We perform a 2D Monte Carlo simulation on a square 
lattice. (Equivalently, the simulation can be described 
as an SOS model in which the height of the surface is 
either or 1.) The energy barrier for a move is taken 
to be proportional to the number of in-plane nearest (n) 
and next-nearest (n') neighbors of the atom before hop- 
ping. Next-nearest neighbors are included to reduce the 
anisotropy of the steps. Moves leading to an adatom 
(i.e. an atom with no occupied nearest or next-nearest 
neighbor sites) are forbidden. Atoms are allowed to hop 
only to nearest-neighbor sites. Electromigration is taken 
into account as a direction-dependent bias in the hop- 
ping barriers. The total barrier energy is then written: 

Eijk = Sdiritj +71^,) +ecmCOs([fc - l]7r/2 - ^o), (32) 



where i and j are position coordinates along the x and y 
axis, respectively; Uij and n'^j are the number of nearest 
and next-nearest neighbors, respectively; k = 1,2,3,4 is 
the direction of the move (to one of the nearest neigh- 
bors); and Ed and Scm are the energies associated with 
diffusion and electromigration, respectively. Explicitly, 
£cm — -Fsta/2; 6*0 is the angle between of electromigra- 
tion force and the [10] axis of the lattice. Atoms and 
moves are picked randomly, and a move is performed 
with probability proportional to exp(— iJij^/fc^T). 

A circular vacancy is chosen as the initial state. The 
late stage of the instability has been the subject of a re- 



cent controversyrJiia As for insulating voids in Ref. |2^, 
our atom and vacancy islands (equivalent to conduct- 
ing voids) split when they are unstable, as illustrated in 
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Fig. 3b. For weaker electromigration (i.e. smaller 
the cluster is stable, as shown in Fig. 3a. In this figure, 
the steady state of the cluster at long time is slightly 
elongated along the electromigration axis. This elonga- 
tion is probably a consequence of the anisotropy of Dl 
or /3, which is not completely avoided in the simulations. 

Quantitatively the instability appears for radii larger 
than Rc ~ 13a at fcsT = 0.6£, where the characteristic 
energy e is the energy of a single bond. We note that an 
elementary kink costs energy e in our model. In the "re- 
stricted" approximation, in which only elementary kinks 
(i.e. those repositioning the step by a single row) are 
allowed, the step stiffness becomes: 



13 = {2 + e^/^^'^) 



ksT 
2a 



(33) 



Using this relation, we find fi = 2.2e/a. Now from the 
criterion |x| > Xc we get Rq > Rc ~ 10.8a, in good 
agreement with our simulations. Note that the insta- 
bility can also be used to determine (3 or F when the 
other parameters are known. From our simulation, with 
Rc = 13a, we find = 3.3£/a. 

As expected from the model, the same instability is 
found for atom islands. 




• 



FIG. 3. Monte Carlo simulations of a vacancy island in the 
PD limit, with F = O.le/a and fc^T = 0.6e. a) Stable island, 
at Ro = 10a. b) Splitting, at Ro = 15a. 



IV. 2D TRANSPORT REGIME 

A. Steady states 

In the opposite limit of 2D transport, the equilibrium 
behavior of these islands is well knownJlj The cluster 
diffusion constant is 



7ri?o V -Ro 



1 



1 



Ro 



(34) 



The corresponding scaling expectation of Eq. (^) is con- 
firmed readily for weak electromigration, where we can 
work to first order in The islands remain circular 

and drift with velocity 



V 



Ro 



Ro 



Ro 



Ro 



(35) 



The first conclusion is that an Ehrlich-Schwoebel effect 
((i_ > d+) is needed for non-zero drift. We note in par- 
ticular the remarkable prediction that for EC, the drift 
velocity actually increases with increasing island size. 
This behavior has been seen in Monte Carlo simulations 
by Wickham and Sethna,Q and leads to an exponential 
coarsening of an assembly of clusters when coalescence is 
induced by drift- velocity differences. Eg. (^) allows us 
to evaluate d+ inside the voids of Ref ^. Their findings 
lead to the expected result that the length d+ increases 
as temperatures decreases. Moreover, their phenomeno- 
logical fitting expression for asymptotically large islands 
(cf. the caption of their Fig. 6), \V\ = uq — wi(a/i?o), can 
be retrieved from an expansion of Eq. (|3^) for Ro ^ d+ 
(and d- — > oo). Thus, we can identify their coefficients: 
vo = a^DCgg/fand vi = vod+/a. 

Since Eq. (|3^) does not depend on cf), both atom and 
vacancy islands drift in the same direction in the 2D 
transport regime. Thus, we obtain a very simple way to 
identify the dominant mass transport mechanism, since 
atom and vacancy islands drift in opposite directions in 
the PD limit. (Cf. Eq. (|l|).) Note also that when d- > 
d+ — expected for the typical Ehrlich-Schwoebel effect 
— the direction of motion of the (adatom or vacancy) 
cluster in the 2D transport regime is opposite to the 
electromigration force. (Cf. Fig. 2.) 

For a vacancy or an atom cluster in the TD limit, the 
circular steady state is an exact solution of the equations 
of motion (even in the nonquasistatic limit). As soon as 
the attachment is not instantaneous {d+ ^ 0), this is 
no longer true. To second order in a non-circular 
shape is found in steady state for the "interior" model 
(e.g. appropriate to a vacancy island with an infinite 
Ehrlich-Schwoebel barrier) when d+ ^ 0. This shape 
is elongated perpendicular to the electromigration force 
for vacancy islands (along the electromigration force for 
atom islands), as can be understood from an intuitive 
argument: For small sticking probability, atoms detach- 
ing from the sides of the vacancy island, will drift along 
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the electromigration axis, but will also have a residual 
drift toward the center of the vacancy if they attach only 
after several attempts. 

A macroscopic description, analogous to the Bernoulli 
effect, provides an intuitive understanding of this phe- 
nomenon. According to this effect, the pressure variation 
perpendicularly to a flow is proportional to the kinetic 
energy involved in this flow. Considering the flow of 
adatoms inside the island, this relation takes the form: 



1 

Jatt X mVatt oc -mc, 



eq ^ adatJ 



(36) 



where the l.h.s. is the pressure variation and the r.h.s. 
is the kinetic energy of the flow. The adatom mass di- 
vides out, it is irrelevant here because the dynamics is 
overdamped. The first factor on the l.h.s.. 



Ja 



(37) 



is the flux of atoms attaching to the steps. The second 
factor mVatt is the momentum associated with atoms 
attaching to the step, i.e., the product of the adatom 
mass m and the macroscopic velocity 



D/{Ro 



(38) 



associated with motion of an adatom across the island 
during its lifetime. On the r.h.s. of Eq. (|36|), we use 
Vadat, defined as the adatom drift velocity needed to 
move the step at velocity V. Hence, 



VadatCeq = V / a'^ 



(39) 



where the cluster velocity V is given by Eq. (^5|) As a 
macroscopic property, Vadat is different from the drift 
velocity (v) of one adatom as defined in Eq (|l|) . Using 
Eqs. (p^-B9) in the Bernoulli relation Eq. (|36|), we get 
the effective concentration change induced by the elec- 
tromigration force: 



i?0 rf-t 



(40) 



Moreover, the concentration variation associated with a 
change of local radius is known from Gibbs-Thomson 
relation 



C Cgg CX Cgg 



TSR/Rl 



(41) 



Balancing the destabilizing effect in Eq. ( pO| ) to the sta- 
bilizing one in Eq. (|4l|), we get a condition for a steady- 
state: 



5R 1 -Rq d-^-Ri) 



Ro T e d++ Ro 



(42) 



Expanding our interior model to second order in 1/^, 
the island deformation involves only the n = 2 mode. 
Defining 



A ^ :r^[p(V2) + p{-n/2) - p{0) - p{n)] 

ZHq 

we find, quantitatively, 

1 Rq Rq 



A = 



r 6^2 + Ro 



(43) 



(44) 



in agreement with the evaluation in Eq. ( p^ ) based on 
the Bernoulli effect. In the EC limit, with Rq ^ lO'^a 
and ^ ~ lO^a, one finds that A ~ IQ-^. Thus, the 
deformation is small for an electromigrating island with 
interior dynamics. 

When mass exchange with the exterior is allowed, 
there is no steady state. To investigate further, we study 
the case of an (atom) island exchanging matter only 
with the exterior, but with desorption allowed. A term 
Feq — c/ts is added to Eq. (^). Although the island is 
not stable in this case, we further consider its behavior in 
order to analyze the tendency to deviate from circular- 
ity as it shrinks. To first order in 1/^ we find a circular 
steady state with velocity 

where Kn = Kn{Rvil Xg) is the modified Bessel K func- 
tion and Xs — {DtsY^^. The n = mode is unstable, 
and the island has a characteristic decay time 



Tde 



AttRI XsKo + d+Ki 



(46) 



To second order in ^ ^ the island shape is not circular 
anymore, and the resulting deformation reads: 



A 



d+RQXs 2xs Rq{K2-KiKz)-XsK2Ki 
'Sr 8^2 Ki^K:i XsKi+d+{KQ + K2)/2 



(47) 



Since A is negative, the elongation of an atom island is 
now along the electromigration axis. 

This can also be qualitatively understood from the 
Bernoulli effect, and the origin of elongation along the 
electromigration axis can be traced back to a change 
of sign of the curvature in the Gibbs-Thomson relation 
(pl|). To adapt the local equilibrium at the step to the 
lowering of the pressure on the sides of the island, the 
curvature has to be locally decreased, so that the island 
elongates along the electromigration axis. In the limit of 
small desorption, we get: 



A = -- 



1 



2r e Ro 



(48) 



Note that A diverges when oo, in agreement with 

the previous claim that there is no steady state in this 
limit. We notice that this formula is similar to Eq. (jij), 
with a factor of Rq replaced by the new cut-off length Xg ■ 
Taking r - a, d+ ~ W^a, Rq ~ lO^a, x, ~ lO^a, and 
^ ^ lO^a, we get A ~ 0.5. Hence, a large deformation of 
the island can be found in the exterior model. 
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B. Validity of the quasistatic limit 

In this section, we consider a vacancy island in the 
absence of electromigration and in the non-quasistatic 
regime. We use polar coordinates r and 9. The nth 
component of the Fourier transform of the concentration 
with respect to 9 is expanded as 



(49) 



Using this expression in the interior model, we look for 
solutions of the form anm(t) = exp{iujt)anm- By sym- 
metry a„m = a,-nm, and from analyticity of c at r = 0, 
o-nm = when n — m is odd. Using the recursion equa- 
tions resulting from Eqs.(0j|), it is the possible to show 



that ar, 



when m < Inl, and that all a„m can be 



calculated from the free parameters a„[„|. We thus cal- 
culate the anm as a function of the a„|„|, and plug this 
result into the boundary conditions Eqs. ([7[]lO|) -with 
J St = 0. We find the non-quasistatic dispersion relation 
from the requirement that the prefactors of the a„|„|'s 
should cancel: 



much smaller than that of the adatoms on the terrace 
(r.h.s.). 

The mode-dependent term on the r.h.s. of Eq. ([5^), 
related to attachment-detachment, indicates that the 
quasistatic limit is more difficult to achieve at large 
|rt|, i.e. small wavelengths, because these modes have 
shorter relaxation times. Using a small-wavelength cut- 
off WC ^ Ra/a on the r.h.s. of Eq. (|5^), we see that 
Eq. ( p3| ) is valid in the limit of slow kinetics, <^ D j a. 
Conversely, when Vj^ ^ -D/a, we find 



aVceq < 1 , 



(54) 



which is the criterion for the quasistatic limit to be valid 
in the fast-kinetics (i.e. TD) regime. 

As an example, let us consider the case of Si(lll) at 



T - lOOOK. From Ref. y, we have a 
d+ — D/v+ ~ lO^'a, and from Ref. 
( ^3|) is then easily checked to be valid: 
approximation can be used. 

C. Stability analysis 



-eg 



< 0.1 and 
F - lOA. Eq. 

the quasistatic 



A = -a^Ceq — {n^ - 1) X 



(50) 



where /„ is the modified Bessel / function, and c^q — 
Cg^exp(— F/i?o)- Since A = iojRl/D, 5Re[A] is propor- 
tional to the growth rate of the perturbations. In the 
quasistatic limit, Eq. (^) is solved on the interior terrace 
instead of Eq. (^. Wc then find the following dipersion 
relation: 



A = -a^Ceq — (n^ - 1) 
Ho 



d+ 
Rq 



(51) 



We now go back to the electromigration problem, and 
perform a linear stability analysis of the non-quasistatic 
model equations. We consider the TD limit {d+ = 0) 
of the interior model. We perform a stability analysis 
of perturbations with respect to the steady state with 
constant concentration c = c^q and circular shape R = 
Rq. To do so, we again expand the concentration using 
Eq. (^), now in Eqs. (@-||), and seek solutions of the 
form anrn{t) = exp(ia;i)a„„i. Matching coefficients for 
for each r™ leads to a set of recursion equations for the 

Inrn's for TTi > 2: 

iU! / 2 2\ 1 [/ ^ 

-j^an,m-2 = -anm(n - m ) ^ '^[\^ + m)an+l,m+l 



{n — m)an- 



l,m-l 



(55) 



Note that in the quasistatic limit, the n — mode is 
frozen (i.e. A = when n = 0). Moreover, we obtain 
here the standard TD (or EC) limit by letting ^ 
(or d+ oo), leading to iu) ^ R^^^ (or ~ Rq^)- From 



Eq. (|5C|) , the quasistatic limit Eq. ( |5l| ) is recovered when 
|A| ^ n?. Using the quasistatic expression for A, this 
condition can be rewritten: 



2~ -^0^^+ 
1 a CeqV+ <C I I 



D 



(52) 



where \n\ > 1. For the approximation to be valid for all 
modes, it is sufficient to satisfy the stronger condition 



Fa^Ce 



(53) 



This inequality should be interpreted as follows: the ef- 
fective diffusion constant of the step (l.h.s.) must be 



supplemented by the relations 

Ono = for n 7^ , 
o-ni = for |ri| 7^ 1 . 



(56) 



Boundary conditions Eqs. ([7[]l0|) then provide two set of 
relations for the a„„i's and the p„'s. When the p„'s are 
eliminated, a set of relations among the a„m's is found: 



= 



OO 

m=0 



iujR^ 
D ' 



Rq 
"2C 



— (n^ - l)T\nianni 


(a„+l,m + an-l,m)| • (57) 



Numerical solution of this set of equations reveals an 
instability for Rq/^ > 0.1 related to the existence of a 
positive eigenvalue iui in this range of parameters. It is 
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interesting to note that the instabihty threshold is in- 
dependent of r. To investigate the nature of the insta- 
bihties of these islands, we again perform Monte Carlo 
simulations. Rather than the splitting seen in the PD 
limit, we find a new type of instability. 

D. Monte Carlo Simulations 

The code used here is similar to that used for a vacancy 
island in the PD limit, but with adatoms now allowed in- 
side. Moreover, in order to tune attachment-detachment 
kinetics, the hopping probability is decreased by a fac- 
tor p, < p < lp-.when the number of neighbors changes 
during the hop.LJ 




FIG. 4. Monte Carlo simulations of a vacancy island in the 
TD limit, with F = O.le/a, ksT = 0.6e. a) "Stable" island, 
at -Ro = 10a. b) Slit formation, for Rq = 30a. 

At high temperature fc^T > 0.7e, micro-vacancies es- 
cape noticeably from the steps. This might heal possible 
instabilities by producing a micro- vacancy where a cusp 
is present. We therefore chose to perform our simulations 
at ksT = 0.6e. Even though a few micro- vacancies still 
escape from the clusters at this temperature, they are 
rare enough so that their effect on the dynamics can be 
considered negligible. 

No drastic change in the island shape is seen at the 
threshold given by the linear stability analysis. Let us 
define the front and the back side of the island with re- 
spect to cluster motion as the first and the last part 
passing at a given value of x. When electromigration is 
increased, the front side is stabilized. The back side is 
destabilized and exhibits chaotic behavior, but its rough- 
ness remains finite, as shown in Fig. 4a. We conjecture 



that this is the instability that we found from linear anal- 
ysis. Although some instability is seen in this regime, we 
classify the island as stable because its shape is not af- 
fected drastically. 

Increasing the electromigration force, we see that un- 
stable slit-like shapes appear, as depicted in Fig. 4b. 
Since these shapes are seen for all orientations of the elec- 
tromigration force, they are not due to lattice anisotropy. 
This instability contrasts qualitatively with the splitting 
found in the PD regime (see Fig. 3b). 

Changing the value of p does not seem to change this 
scenario, but we have not performed simulations in the 
limit where p is small {p < 0.1) due to the large computa- 
tion time needed there. Numerical solution of the model 
equations would be needed to understand in greater de- 
tail these different regimes. 

V. ELECTROMIGRATION BIAS IN 
ATTACHMENT-DETACHMENT 

Finally, we introduce a possible electro-bias between 
attachment and detachment. Extracting one atom from 
the step, not counting any bond breaking, requires work 
aFem'fi along the normal to the step.a From a linearized 
Gibbs-Thomson relation, we find that the effective equi- 
librium concentration becomes 

Ceg = c%{l + aF,^-h/kBT) . (58) 

Defining = Fem-n/fesT" on the lower or upper side 
of the boundary step edge, respectively, we find that 
electro-bias contributes to the velocity of Eq. ( p5| ) an 
extra term 

where the subscript indicates the n — I component of the 
Fourier transform of l/£,±{0). When there is no Ehrlich- 
Schwoebel effect {d+ = d-), V from Eq. ( ^5| ) vanishes, 
and the drift velocity is only due to this new contribution 
combined with PD effects. In the EC limit, 6V does not 
depend on the size of the cluster, while in the TD limit 
SV ^ l/i?o- This is the same size dependence as that 
associated with the PD limit. Since SV decreases with 
increasing island size, we expect this contribution to be 
small for large islands. 

The deformation of the cluster is proportional to 1/^, 
rather than 1/^'^ as seen in the other regimes. Thus elec- 
trobias might induce larger deformations of the clusters. 
However, since we do not know the angular dependence 
of these forces, we cannot go further in the analysis. 

VI. SUMMARY AND PERSPECTIVES 

We have studied the response of single-layer clusters 
to electromigration in different regimes associated with 
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different mass-transport mechanisms. Some of the main 
resuhs are collected in table 1. The drift velocity and 
shape changes agree with those found in earlier published 
calculations ,□ and give a theoretical explanation and ex- 
plicit expressions for those found with Monte Carlo sim- 
ulations. 

The instability is found to depend qualititatively on 
the mass transport mechanism. PD is associated with 
splitting, and 2D transport with slit formation. Al- 
though the instability threshold of the latter may be be- 
yond the experimental range for 2D islands, it could well 
be important in the case of disordered |Cix}igen domains in 
YBaCuO,crQ and for void dynamics J3'E3K3 where higher 
current densities are used. In the latter case, 2D trans- 
port might occur on the substrate across the void. 

Fluctuations are found to be affected by electromigra- 
tion. The non-equilibrium cluster diffusion constant and 
roughness were studied. The island roughness is found 
to diverge as the instability threshold is approached. 

Non-linear analysis is needed for a better understand- 
ing of these instabilities. Moreover, cross-over regimes 
should be studied in order to catalogue completely the 
behavior of these islands. For a more realistic description 
of the islands, especially aJ. l o w . temperature, anisotropy 
should also be included.a'E§Ea Indeed, it was already 
shown experimentally that extremely anisotropic prop- 
erties (such as diffusion) could lead to very different 
behavior .E3 





PD 


TD 


EC (AD) 


Length 
criterion 


Rl < 4(i?o + d) 


-Ro 3> d+ , d- 


Ro < d+ , d- 


{V{Ro)} 


1/Ra 


1 


Ro 


Atom / vacancy 
direction 


Opposite 


Same 


Steady state 
shape 


Circular 


None 
Non-circ. 


Instability 
morphology 


Splitting 


"Slit" II F 


Slit ± F 






Anisotropic 



TABLE I. Summary of results. In general the entries ap- 
ply to both atom and vacancy islands. The criteria for the 
three regimes in terms of the characteristic lengths are given 
in the first row. The inequality for PD is dominant; when it 
holds, one has PD even if the inequality for TD or EC is satis- 
fied. The characteristic length £c is Dl / aDc^q, with equilib- 
rium values for the parameters. For the third row, the relative 
direction of the drift velocity of atom and vacancy islands is 
tabulated. In the fourth and fifth rows, the upper (lower) 
entry is for atom/exterior (vacancy/interior) islands. PD, 
TD, and EC indicate the mode of mass transport: periphery 
difi'usion, terrace diffusion, and evaporation-condensation, re- 
spectively. Reidentifying the letter D, the last is often called 
AD (attachment-detachment) to indicate that the process in- 
volves just a 2D "gas" on the terraces. 



This study is a step toward the understanding of the 
evolution of more complex surface geometries in the pres- 
ence of a diffusion drift, such as 3D islands or groves or 
assemblies of clusters. To complete this exploration, the 
fluctuations in the 2D transport regime should be stud- 
ied. It would also be useful to check with Monte Carlo 
simulations the predictions about the diffusion constant 
and the divergence of morphological fluctuations. 

We note in closing that, as the surface is affected by 
electromigration, so electromigration itself is affected by 
the morphology of the surface. Step shadowing, adatom 
crowding, or surface roughness induce changes in the lo- 
cal electromigration force.&ll A quantitative description 
of experiments should take these feedback mechanisms 
into account. 
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